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Abstract 

The definitions of the main notions related to the quantum inverse 
scattering methods are given. The Yang-Baxter equation and reflection 
equations are derived as consistency conditions for the factorizable scat- 
tering on the whole line and on the half-line using the Zamolodchikov- 
Faddeev algebra. Due to the vertex-IRF model correspondence the face 
model analogue of the ZF-algebra and the IRF reflection equation are 
written down as well as the ^-graded and colored algebra forms of the 
YBE and RE. 
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1 Introduction. 



Integrable and/or solvable models have been always quite important for the- 
oretical physics. One of their attractive features is the direct possibility to 
go beyond perturbation theory giving a solid background for theoretical hy- 
pothesis and constructions. Another valuable characteristic is related to the 
combination of different mathematical methods required to solve some specific 
model thus giving rise to mutual interrelations among formally separated fields 
of mathematics. The start of the recent activity in the quantum integrable 
models (see reviews [I], |2|, ||, [|, [|, ||, [7|, |8|] and Refs. therein) was definitely re- 
lated with the development of the soliton theory, although the influence of the 
preceding pioneering contributions is out of discussion (see [[TO], [II], [T2|, |T3J and 
Refs. therein). This is a large field of research and even recent monographs 
TT], [TJ, [TJ, [TJ, |TJ, [T§ [T7|, |TJ do not overlap much in between. 

There is variety of quantum integrable models and quite a few interrela- 
tions among them. In these lectures we restrict ourselves to those models 
the integrability of which is related more or less directly to the Yang-Baxter 
equation and to reflection equations. More of that, we will discuss mostly 
the algebraic properties of these equations, their solutions and corresponding 
integrable models putting aside very elaborated analytical techniques and/or 
problems (of great physical significance) such as the thermodynamic limits of 
finite size systems, massless and/or conformal field theory limits, continuous 
limits of the lattice models, correlation functions and form factors, critical ex- 
ponents etc. The main notions related to the Yang-Baxter equation (YBE) will 
be introduced: factorizable scattering on the line, the Zamolodchikov-Faddeev 
algebra, the fusion procedure or the bound state scattering, the Yang-Baxter 
or transition matrix algebra, integrals of motion. The factorizable scattering 
on a half-line gives rise to the reflection equation (RE) and a boundary oper- 
ator. Most of the applications of the YBE can be extended to the RE case 
with appropriate modifications. Due to the vertex-face model correspondence 
of statistical mechanics we introduce an interaction round face analogue of the 
Zamolodchikov-Faddeev algebra and the corresponding reflection equation for 
the boundary weights. The super or Z 2 -graded analogue of these constructions 
with the generalizations to the color algebras will also be given as well as some 
particular integrable systems with finite degrees of freedom. The interrelations 
among few forms of the Bethe Ansatz: coordinate, algebraic, analytical and 
functional will be mentioned as well. 

More particular applications of the Bethe Ansatz technique to specific inte- 
grable models with direct physical meaning ( a high energy limit of the QCD, 
the Azbel - Hofstadter problem ) and of the reflection equation formalism to 
the Chern-Simons theory and quantized moduli spaces of flat connections can 
be found in other lectures of this School. 
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2 Yang- Baxter equation and Zamolodchikov- 
Faddeev algebra 

Let us start the presentation which will be full of definitions and technicalities 
from a specific integrable model to have in mind an example. This is the 
one-dimensional Bose gas consisting of n sort particles with the Dirac delta- 
function two particle potential. The field operators ip a (x),ipl(y) satisfy the 
canonical commutation relations, a,b — 1, 2, • • • , n 

[4> a (x),4(y)] = $a b 5(x-y) , [Mx)My)] = o- (!) 

The field Hamiltonian is 

H= / dx(d x tpl(x)d x ip a (x) + cipl(x)ipl(x)ip b (x)ip a (x)) (2) 



The translational invariance (on the whole line or with the periodic boundary 
conditions, circle) and the internal U(n) invariance are obvious. The corre- 
sponding symmetry generators are 

P = -i J dx^l(x)d x ^ a (x), Uab = J dxi>l(x){u ab ) cd i> d (x). (3) 

where u a b are n 2 generators of the unitary Lie algebra u{n). One of these 
operators is the particle number operator 



N 



f dx^ a {x)^ a {x) . (4) 



The Fock space of states is the direct sum of the number operator iV 
eigenspaces 

oo 

U F = n M (5) 

M=0 

The common eigenfunctions of operators N, P, and H (2) 

. M 

q, M = d M x^(l,---,M\X 1 ,.--,X M )H^ ai (x i )\0) (6) 

i=l 

are constructed as appropriate linear combinations of the one particle eigen- 
states (the plane waves). The numbers j = 1, 2, • • • M in the argument of ^ 
(6) refer to both coordinate Xj and isotopic aj indices. The coefficients in the 
linear combination depend on the one particle parameters (momenta) Xj and 
on the isotopic indices aj 

M 

= E MM> {K})ex P (i J2 Km X m ) (7) 

where a are all elements of the permutation group Sm- 
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Such a form of the Hamiltonian eigenfunctions is known as the coordi- 
nate Bethe Ansatz. The conditions of the wave function continuity and its 
appropriate derivative jumps on the hyperplanes Xj = Xj + i (the sewing condi- 
tions) define the coefficients A CT ({aj}|{Afc}). The coefficients A a and A a / with 
a' = o-jcr where <jj is the transposition of the indices j,j + 1, are related by 
the two particle ^-matrix [20]: 

A a , = S(Xj - X j+1 )A a , (8) 
where for the model chosen (2) 

S{\) = (A + icP)/(X - ic) , (9) 

and V is the permutation operator in C n Cg> C n . The periodicity condition for 
the system on finite interval (0, L) results in the Bethe equations for the set 
of M momenta Xj 

exp(iXjL) = - ]J s jk(*j - A*) , (10) 
k 

where in the ordered product k — j + 1, M — 1, M, 1, — 1. Hence the 
meaning of the RHS is the scattering matrix of the j-th particle on the other 
(M — 1) particles. This would be just a phase factor for the scalar particles 
(S(X) = (A + ic) j (A — ic)) , but for the n component case one has to diagonalize 
the complicated scattering matrix to arrive to a system of scalar equations. All 
quantities in the RHS are particular values at A = Xj of the transfer matrix 
(fc = l,2,..,M) 

t(A; {A m }) = tr a T(X; {X m }) = tr a Y[S ak (X - X k ) (11) 

k 

of the inhomogeneous GL(n)-spin magnet of M sites with the transition or 
monodromy matrix T(A; {A m }). The trace in the expression for the transfer 
matrix t(X; {A m }) is taken over the auxiliary space V a = C n , while t(X; {A m }) 
is an operator (matrix) in the space Ilfcli(C n )fc- The important property of 
the commutativity of the transfer matrix for different values of the spectral 
parameter 

[t(X;{X m }),t(fi;{X m })] = 

follows easily from the fundamental commutation relation for the transition 
matrix T(A; {A m }) (see (15)). 

Hence the RHS of (10) t(Xj] {A m }) can be diagonalize simultaneously. As 
result the hierarchy of the Bethe Ansatze appears (or the nested Bethe 
Ansatz) and the complete parametrization of ^({jjllAj}) has n sets of "quasi- 
momenta" including Xj (see below). 

The consistency condition of this system is the Yang-Baxter equation 
(YBE) for the S-matrix S{X): 
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Sjki^j — ^k)Sji(\j — Xi)Ski{Xk — X{) — 



— SkiyXk — Xi)Sji(Xj — Xi)Sjk(Xj — Afc). (12) 

One uses very often in the general setting the notation R ( the i?-matrix ) 
for the YBE solution. 

The complete scattering matrix 5'({Afe}) of the M particle is given by the 
ratio of the coefficients of the incoming wave (Ai < ... < Am) and the outgoing 
wave 

exp(i^2XjXj), exp(i^2X M -j+iXj) ■ 
j j 

This ratio is factorized into the ordered product of M(M — l)/2 two particle 
S'-matrices (9) ( the factorizable scattering on the line ). 

It was proposed for an algebraic description of the factorizable scattering 
in the general case to introduce a set of (annihilation) operators Z a (X) fll9t 
satisfying the commutation relations ( the Zamolodchikov algebra ) 

Z a (X)Z b (u) = S ab , cd (X - v)Z d (u)Z c (X), (13) 

where S(X — v) is an n 2 x n 2 matrix. Using the associativity property of this 
algebra and changing the order of the product Z ai (Xi)Z a2 (X2)Z a3 (X 3 ) in two 
possible ways one arrives to the consistency condition (12). Extending the 
Zamolidchikov algebra (13) by adding n more (creation) conjugated operators 
Z\[v) one gets the Zamolodchikov- Faddeev algebra (ZF-algebra) 

Z a (X)Z\{v) = 5 ab 5(X + Z\{v)S acM {X - u)Z d (X). (14) 

It is useful to write down the ZF-algebra in a compact matrix form by 
introducing the n component column A(X) = (Zi(A), Z n (X)Y and the n 
component row A\v) = (Z\(u), Z^(y)). Then the defining relations of the 
ZF-algebra are 

A(X) ® A(u) = A 1 (X)A 2 (u) = S 12 (X - u)A 2 (u)A 1 (X), 
A\(X)aUu) = A\{v)A\{X)Sl{p ~ A), 

A 1 (X) ® A\(u) = h6(X -u) + A\(v)S n {v - A)A 2 (A), 

where subscripts refer to the corresponding isotopic spaces C n ® C n = ViV 2 
and 5*21 = VSi 2 V, S 12 = VS\ 2 . Due to the unitarity property of the S- 
matrix: 5i2(A — v)S 2 \{y — A) = I\ 2 and the YBE one can construct the Fock 
space representation Tip of the ZF-algebra using the generalizing symmetrizing 
operators which include S instead of the permutation operators [37]. 
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To solve the non-linear Schrodinger equation (NS) ( the Heisenberg equa- 
tion of motion for ip a (x,t) with the Hamiltonian (2) ) in the framework of 
the quantum inverse scattering method (QISM) [1 - 7] one has to find the 
corresponding auxiliary linear problem 

-^-T(X,x) = L(X,x)T(X,x) . 
ax 

For the model in question L depends on the spectral parameter A and the 
original dynamical local variables (1) 

L(A, x) = XJ + ^(ctfj a (x)e at n + i + ipl(x)e n+1 , a ) 

a 

where J is (n + 1) x (n + 1) diagonal matrix: J = diag (I n , —1) and are 
the (n + 1) x (n + 1) basis matrices. Hence the auxiliary linear space here 
is C n+1 . This is the L-operator of the QISM or the classical soliton theory 
[13]. Solution to the operator valued matrix first order equation with normal 
ordering with respect to the local fields (1) and the vacuum |0) 

T(A, x) =: exp( f L(X, y)dy) : 




defines the transition or monodromy matrix T(X,x). Its entries are the 
new variables of the model (the quantum scattering data (QSD)). 
For the NS case it is natural to represent T(A, x) in the block form 

T = 

where A is an n x n matrix, B and C are n component vectors, D is a scalar. 
The quadratic commutation relations of the new variables are defined by the 
fundamental relation |], [| 

R(X - v)T x {X)T 2 {v) = T- 2 (u)T 1 (X)R(X - v) (15) 

where the standard QISM notations [1-7] are used T\ — T eg) /, T 2 = I <8> T to 
embed matrices in C n into C n <8> C n . From the structure of the L-operator one 
can conclude that the n elements of the row C(A) act on the vacuum as some 
creation operators. This is the starting point to construct the eigenf unctions 
of the transfer matrix t(^) = trT(^i) = trA + D algebraically: 

M 

^/ = n^( A J )°) 

3=1 

in the framework of the algebraic Bethe Ansatz [|I|, 0, H]. 

The mentioned above commutativity of the transfer matrices t(^) = trT(fi) 
follows from (15) taking the trace of T{T 2 = Ri 2 1 T 2 TiR 12 over both spaces 
C n ®C n . The trace form of the integrals of motion generating function t(n) for 
continuous as NS (2) or chain models as (11) gives rise to the periodic boundary 
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conditions. For the M site spin Hamiltonian H = (Snll 1 ^n,n+i) + ^m,i the 
M-th site spin interacts with the first one. The treatment of the non-periodic 
boundary conditions in the algebraic framework of the QISM requires the RE 
[29, 32, 34, 31] (for the coordinate Bethe Ansatz see [12, 28]). 

Using block decomposition for the i?-matrix one can rewrite the compact 
form of the fundamental relations (15) in terms of the blocks A, B, C, and D. 
Restricting to the case of finite number of particle and after appropriate limit 
to the whole line one gets the ZF-algebra realization in terms of the QSD [50] 

Z(X) = (A)-\X)B(X), Z\v) = C{y)D-\v) . 
The i?-matrix of (15) satisfies the YBE 



R 12 (\ - aO-Ri 3 (A - v)R2a{}i -v)= R 23 (n - v)Ria{>< - v)Ri2^ - fj) ■ 

For the NS the i?-matrix has the same structure as (9) (but it is (n + l) 2 x 
(n + l) 2 matrix). 

In the general situation the i?-matrix depends on the spectral parameter u 
and some other parameters R(u;r), ...). Although there is no complete math- 
ematical theory of the Yang-Baxter equation, variety of solutions are known 
as well as different fields of their applications. In particular, many solutions 
are related to the simple Lie (super) algebras. They are classified by the Lie 
algebra, its irreducible representations Aj and the spectral parameter depen- 
dence: rational, trigonometric and elliptic ones [3, 9, 11, 48]. The sl(2) spin 
1/2 R- matrix related to the AXA-magnet is used in few lectures of these Pro- 
ceedings. The recent development relates the spectral parameter dependent 
i?-matrices with the affine Lie (super-) algebras fllEf . There are also solutions to 



the YBE with the spectral parameter on the algebraic curves of higher genus ( 
the Potts models ). The i?-matrices acting in infinite dimensional spaces can 
be found in papers [58]. 

Some particular properties of the YBE solution R(u; f], ...), which are im- 
portant for different applications ( but not always valid for a given solution ) 
are: regularity 

R(0) = <P( V )V , 

P-symmetry 



T-symmetry 



unitarity 



VR 12 (u)V = R 21 (u) = Ru(u) 



R T 12 {u) = R 12 {u) 



R 12 {u)R 2 i{-u) = p{u)I 
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crossing symmetry 



R\ 1 2 (u)R%(-u-r ] ) = au)I, 



R(u; rj) = I + r)r(u) + 0{tf 



quasiclassical property 



where r{u) is the classical r-matrix [g, [13j and 4>(u), p(u), £(u) are some 
functions related to the i?-matrix normalization. Many i?-matrices have only 
PT-symmetry: Rj 2 {u) = i?2i(u). The regularity is used to extract from t{u) 
of the lattice models the integrals of motion which are local in terms of initial 
spin variables when the i?-matrix itself is the L-operator. The quasiclassi- 
cal property gives rise to the direct connection of the quantum model to the 
corresponding classical one [|3], [TBI . 

It is easy to see that the product of the .R-matrices Rn{ u )R2^{v) an d 
Rvziw) are intertwined (like (15)) by the matrix Ru/(u — w)R2v(v — w). If 
the original R- matrix Ri2(x) is degenerated into a projector at x = r\ then 
one can project the above product to the corresponding subspace for the fixed 
difference of the spectral parameters u — v = r\. For the Yang solution (9) 
at A = ±zc one has the ( anti )symmetrizer P + , P_ . This is the fusion 
procedure [H] to get new i?-matrices from the known ones. It has a direct 



physical interpretation as construction of the bound state ^-matrix [19], 
Using the fusion procedure the L-operators in the higher dimensional irre- 
ducible representations can be obtained giving rise to the ZF operators for the 
bound states and to the integrable lattice models of higher spins such as the 
spin s XXZ-model. The connection of the i?-matrices and integrable models 
with the simple Lie (super-) algebras is reflected in the structure of the Bethe 
equations: they include r sets of " quasimomenta" , where r is the Lie algebra 
rank, and the Cartan matrix [37, 47]. 

Omitting the spectral parameter dependence in the YBE one gets still very 
interesting equation solutions of which ( the constant i?-matrices ) can be 
considered as structure constants of the quantum groups and the quantum 
algebras. Then the quadratic relations (15) ( without the spectral parameter 
) are the defining relations on the n 2 generators T a ^ of the corresponding 
quantum group. The constant i?-matrices Rij have also direct relation to 
the braid group (BG), for one of the defining relations of the BG generators 
G{ (jj+i Oi = Oi+j <Ji o"j + i coincides with the YBE for R^ = V)\R,)\: 



Ri-li R-ii+l Ri-li — Rii+1 Ri-li R, 



3 Reflection equations and their covariance 

Let us consider factorizable scattering of particles with internal degrees of 
freedom on a half-line [22]. Then even one-particle process in nontrivial (a 
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reflection from the wall) and it is described by an n x n matrix K{u) the 
reflection matrix or the boundary S-matrix. 

For an algebraic description it is useful to add to the ZF-algebra a formal 
boundary operator B with relation [24] 

Z a (u)B = K ab {u)Z b {-u)B . 

Then the two particle factorizability gives rise additionally to the YBE for the 
two body ^-matrix (12) to the reflection equation (RE) 



S 12 {u - v)K 1 {u)S 2l {u + v)K 2 {v) = 

K 2 (v)Si 2 (u + v )K 1 (u)S 2 i(u - v) . 

(16) 

There are natural properties of the reflection matrix K(u) as in the preceding 
Sec. of R(u;r),...) : K(0) = I (regularity); K{u)K(—u) = I (unitarity); 
T-symmetry K T (u) = K(u); the crossing symmetry is more elaborated and 
it involves the S- matrix itself [24]. The constructions of the quantum group 
invariant spin systems [30] uses the RE [31, 32, 42, 54]. Recent field theoretical 
applications of the RE can be found in [23-27, 36]. 

It is interesting to point out that as a consequence of the RE flilf ) the 
quadratic combination 



<P( U ) = ^Z a {-u)K ah Z h {u) 



a.b 



is a "local" field : [<£>(«)> V 9 ( u ')] = EH HI- This property was used to 



construct the metric tensor field for the quantum Liouville theory |52| with 
a slightly different RE and to propose a form of the boundary operator B = 
exp(f£° ip(u)du) in the dual Hamiltonian picture of the boundary conformal 
field theory 0. 

Although for many i?-matrices solutions to the reflection equations were 
found (cf [22, 25, 32, 33] ) there is no direct relation of them to the Lie 
algebra theory. In particular, many of them do not depend on the quasiclassical 
parameter rj (in the fundamental representation for the i?-matrix). 

Let us give two examples of the RE solutions. Due to the GL(n) symmetry 
of the Yang solution (9): [R, M®M], M 6 GL(n) the corresponding i^-matrix 
can be transformed K — > K' = MKM~ X with arbitrary M and the solution 
to the RE is 



K{u) =£1 + uC, C 2 = I. 

In the elliptic case (the eight-vertex R-matrix [11]) the solution K(u) is 
[22, 25, 33] 

/ x(u) y(u) \ 
W I z{u) w{u) ) ' 
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x(u) = 



sn(£ + u) 



y(n) = /j,sn2u(X(l — ksn 2 u) + 1 + ksn 2 u)(l — k 2 sn 2 £ > sn 2 u)~ l , 
w(u) = sn(^ — u) , 

z{n) = fisn2u(X — 1 — (A + l)ksn 2 u)(l — k 2 sn 2 ^sn 2 uY l , 

where snu = sn(u; k) is the Jacobi elliptic function of modulus < k < 1, 
£, A, n are parameters. 

The RE (16) has an important covariance property: if T(u) and K(u) 
satisfy the relations (15), (16) then K'{u) = T(u)K(u)T(—u)~ l is also the RE 
solution provided that the entries of K(u) and T(u) commute [K ab (u),T cd (u)] = 
0. The proof follows easily by the substitution of K'{u) into the RE and using 
few times different forms of the fundamental relation (15) e.g. 

T-\-v)R 12 (u + w)Ti(u) = T^Ruiu + v)T 2 \-v). 
This property gives rise to the Sklyanin monodromy matrix []29], 34 1 



T(u) = T{u)K{u)T(-u)-\ (17) 

If the matrix T{u) is constructed as an ordered product of N independent 
L-operators then T(n) can be interpreted as the monodromy matrix of N site 
lattice model with a boundary condition described by the matrix K{u) or a 
boundary interaction described by the operator valued entries of the matrix 
K(u). To extract the corresponding Hamiltonian and other integrals of motion 
the transfer matrix is constructed using a special trace 

t(u) = trK + {u)T{u)K{u)T{-u)- 1 = trK + {u)T{u) . (18) 

An extra K-matrix K + {u) is any solution of a "conjugated" RE [29, 32, 34, 
51] defined in such way to guarantee the commutativity [t(u), t(w)] = 0. In 
the regular case 

Ran{0) ~ V an , K(0) ~ I , h n , n+1 = {d/duR n , n+1 ) (0) , 

the Hamiltonian is [29] 

M-\ 

H=Y. Kn+i + K[(0) + {tr o K + {0)h M ,o)/trK + {0) . (19) 

n=l 

The structure of the fusion procedure for the reflection equation 

is similar to the i?-matrix case, but the projected combination ( Ruiv) ~ 
projector ) 

K 1 (u)R 21 (2u-r ] )K 2 (u- V ) 
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includes the additional i?-matrix in between the .ff -matrices [32, 51]. 

More general RE with four different i?-matrices ( with or without the 
spectral parameter ) 

R { S K x R { $ K 2 = K 2 K x R$ . 

which are related among themselves by some consistency conditions similar to 
(12) and (15), can be found in different papers these days (see [55, 51] and Refs 
therein ). In particular, this kind of RE was used for the quantum group in- 
variant spin models with topological interaction [42, 43]. Some generalizations 
of the RE are related to the Coxeter groups [44]. 

The constant RE is also of interest as it was in the YBE case. The corre- 
sponding RE-algebras with the if-matrix entries considered as the generators, 
are related to the quantum group homogeneous spaces while the c-number 
solutions K can be considered as representations of a special generator of the 
BG^ in the solid handlebody |yj with the defining relation 



R 12 K 1 R 12 K 1 =K 1 R 12 K 1 R 
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In solvable models of statistical mechanics the i?-matrix defines the Boltz- 
mann weights R a p-^s{u\ rj, ...) of the vertex models, where a, j3, 7 and 5 are spin 
variables on the four edges round a vertex. The models on the dual lattice are 
known as the interaction round face models (IRF). The corresponding Boltz- 
mann weights w(a, b, c, d\u) satisfy the star-triangular equation (relation) 
(STR) or Baxter relation [[Tj| 

w(a, b, g, f\u)w(f, g, d, e\u + v)w(g, 6, c, d\v) = 

9 

= ^2w(f,a,g,e\v)w(a,b, c, g\u+ v)w(g,c,d,e\u) 
9 

The four sites surrounding a face are ordered anticklockwise from the south- 
west corner and u is a complex (spectral) parameter. 

Let us introduce a face analogue of the ZF-algebra: Its generators are 



parametrized by two indices Z a b(u) and satisfy the relations (see e.g. f57j) 
Z ab (u)Z bc (v) = Y w(a, b, c, d\u - v)Z ad (v)Z dc (u) 

d 

with the only one summation over d in the RHS. As in the previous case the 
scattering of a one particle Zd a {u) on N others Zb iCi { v i) w i n gi ye a transition 
matrix 

TV 

T(u,vi, ■■■ ,v N ) = Y[ w{bj,b j+1 ,c j+ll Cj\u - Vj). 
i=i 

by analogy with the YBE and ZF algebra constructions one can formulate 
the following problems: 1) what are the conjugated (annihilation) operators 
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and their commutation relations with Z ab (u)l 2) how to translate the fusion 
procedure to the STR case? 3) what is an IRF analogue to the RE ? 

Let us start from the last question and introduce a formal operator B of 
the boundary which satisfies relations with Z a b(u): 

Z ab {u)B = Q(a, s, b; u)Z as (-u)B , 

where Q(a, s, b; n) is a face analogue of the reflection matrix K. To derive 
an STR analogue of the reflection equation the two ways of transforming the 
product of three operators may be considered 

Z ab (u)Z bc (v)B — > Z ak {-u)Z kd (-v)B 

The resulting relation for the "scattering" matrix w(- ■ ■ ) and the reflection 
matrix Q(- ■ • ) is 

w(a, 6, c, g\u - v)Q(g, c, f\u)w(a, g, f, k\v + u)Q(k } d, f\v) = 



= Q(b, m, c\v)w(a, b, m, s\u + v)Q(s, d, m\u)w(a, s, d, k\u — v) , 

with summation over g, f in the LHS and over m, s in the RHS . 

To relate vertex models with the IRF models the Baxter intertwining vec- 



tors come to play [ ] 1 



Y. R P"-M ( f ) u'(P) z m>i>(v) = J2w(m,m' ,l,l')z m i(X)(j) mm >(a) . 

/3,u rn 

These vectors remind a combination of the operators of the ZF-algebras of 
the YBE and STR types. These vectors can be used to relate the reflection 
matrices of the vertex and IRF models: 

(t>ab{P)Kpv{u)(f>bs{v) = Q{a,s,b;u) , 
where the spectral parameter dependence of the vectors is omitted. 



4 Integrable models with anticommuting vari- 
ables (fermions) 

One of the first integrable models of the quantum field theory-massless Thirring 
model had attracted a lot of attention of theoretical physicists during decades 

H = J dx^di) + g(^) 2 ) . 

2 M.D. Gould informed me that the IRF analogue of the RE was considered by P.A. 
Pearce. 
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It contains Fermi fields ip a {x, t) , ipt(x, t) satisfying anticommutation relations 

{ijj a (x, t),ij>$(y,t)} = 5 aP 5(x - y) . 

The addition of the mass term rmjnp to the Hamiltonian gives rise to the very 
rich dynamical properties of the massive Thirring model. In particular, it is 
dual to the famous sine-Gordon model with the change of the strong coupling 
regime to the weak coupling one. 

The famous (^-potential Bose gas model [20] ( another name is quantum 
nonlinear Schrodinger equation ) (2) was also generalized to the multicompo- 
nent fermionic case in the framework of the coordinate Bethe Ansatz [45] (see 
below). 

Taking into account the " superization" procedure of the theoretical physics 
in the seventies, let us follow this pattern and extend most of the preceding 
equations and constructions to the case of the anticommuting variables. Math- 
ematically it is related with ^-graded vector space, ^-graded algebras or as 
in the theoretical physics text: super-spaces, super-algebras, super-analysis... 

The Z 2 -graded Zamolodchikov algebra is generated by the (m+n) operators 
Zi{u) satisfying 

Zi{v)Zj{v) = (-r^Rij-M (u,v)Z l (v)Z k (u) . (20) 

where i,j = 1, 2, • • • , m + n,p(i) is a parity function: p{i) = (i = 1, • • • , m) 
and p(i) = 1 (i — m + 1, • • • , m + n), hence one has m boson and n fermion 
operators. The coefficients in ( fZUD are c-numbers and the i?-matrix is even: 

p{R im ) = p{%) + p(j) + p(k) + p{l) = . (21) 

Due to the Z 2 valuedness of the function p(i), and the i?-matrix structure 
(^T|) one can change sign factor in (20) to (— )p( fc )pW_ The super- associativity 



requirement for Z a (u)Z b (v)Z c (w) results in the graded Yang-Baxter equation 
(grYBE) 

R ab , al y (it, v)R a , c .j d (it, w)Rv c ,. )hl (V, W) (-)P(«)(PU)+P(a>)) = (22) 

= R bc , Ud (v,w)R ac ,, an (u,w)R a ^ jk ( U) ,)(-)K^(pW+p(«')) , 
where the condition (|2l|) was used to reduce the factors in the LHS and RHS: 

/ _y{a)p{l)+p{a' )p(c)+p(b')p(c') . / _yj{b)p{c)+p{a)p(c')+p{a')p(b') 

Recalling these sign factors one can write the grYBE in the same matrix form 
as previously but now in the Z 2 -graded tensor product of three super-spaces 

R 12 (it, v)R 13 (u, w)R 23 {v, w) = R 23 (v, w)R 13 (u, w)R X2 (u, v) . (23) 
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Considering the third space as an unspecified quantum space (or substituting 
instead of the (m + n) x (m + n) blocks in i? 13 and R23 formal entries) one 
arrives to the graded FRT-relation 

R 12 {u, v)Tx {u)T 2 (v) = T 2 {v)Tx {u)R 12 (u, v) . (24) 



The latter relation has the same sign factors in the component form as (|22|). 
These factors can be interpreted as result of the Z 2 -graded tensor product of 
two even matrices F and G (p(F ac ) = p(a) + p{c)) 

(F <g> G) ab;cd = H Pmj(a)+p(c)) F ac G bd . (25) 

Hence, one gets additional sign factors for T\ and the factor free T 2 (to consider 
more tensor factors the choice of 7\ factor free is more convenient ) 

(TxUcd = (T ® I) ab , cd = (-y^)+^))T ac 5 M , (26) 

(T 2 ) ab;c d — (I ® T) ab -cd = S ac Tbd 

An additional sign factor appears also in the matrix of the permutation (flip) 
operator 

V{v®w) = {-) p{v)p{w) w®v, v,weC mln 

(V) ab ., cd = H p{a)p{b) 5 ad 5 bc . (27) 



To arrive later to the graded reflection equation it is convenient to write ( p6|) 
in a matrix form with the usual tensor product 

(T a ) ar = r(T<g>I)r, (28) 

where {T) ab , cd = (-) p{b)p{a) S ac 5 bd , hence T 2 = I (which is not the case of the 
color algebras). The Z 2 -graded algebraic structures were analyzed also in [35] 
and many R- matrices were found [48]. 

As in the first lecture the simplest integrable models correspond to the 
rational solutions of the grYBE. The GL(m\n) symmetric R- matrix is ( the 
graded Yang solution) 

R{u - v ) = {u - v) + r]P. (29) 
The corresponding models (L-operators) are the GL(m\n) isotropic graded 



magnets and the super-matrix nonlinear Schrodinger equation B^, [53 



L gm (u) = u + riJ2e l3 ® s^(-) p(i) , (30) 
L NS (x, u) =uJ + r] Y^( e ab^ba{x) ± e ba ^ ab (x)) , (31) 
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where Sjj are the generators of the super-algebra gl(m\n), i,j = 1, 2, ■ • ■ , m + 
n; i/j(x), ipl^ix) are bose or Fermi fields according to p{a) + p{b) = or 1 and 
1 < b < N , N+l < a < m+n, J is the block diagonal matrix (1^— I( m + n -N))- 

Among integrable models related to the Z 2 -graded case one can mention 
osp(l\2) - non-linear Schrodinger equation which is described by the ra- 



tional limit of the osp(l\2) i?-matrix and the supersymmetric sine-Gordon 
model with the trigonometric s/(2|l) R- matrix [3]. Although the structure of 
the osp(l\2) i?-matrix is similar to the one of the sl(2) spin 1 i?-matrix the 
solution of the corresponding osp(l|2)-magnet ( or NS ) is obtained [[37) by 
the analytic Bethe Ansatz. Due to the arguments of analyticity, crossing 
symmetry and the bare vacuum the form of the transfer matrix t(u) eigenvalue 
is 

M 

A(u, {vj}) = (u(u - d^HS^u -v k - l/2)5_i(« - v k - 1)- 

k=l 

AI M 

-{{u - \){u - dVHStiu - v k ) - {u{u -d+ l)) N l[S^(u -v k -d), 

k=l k=l 

(32) 

where Si(u) = (u + l/2)/(u — 1/2) and d = 3/2 . The regularity condition of 
A(u, {vj}) in u ( the Manakov's principle ): ResA(u, {vj}) = at u = v k gives 
the Bethe equations for the set of the quasimomenta {vj} 

M 

(5 , i(^fc)) 7V = II s *( v k - Vj)S^i{v k - Vj). 
i=i 

Hence the M super-particle eigenstates are parametrized by one set of the 
quasimomenta according to the Lie super-algebra rank one. 

Another way to define the form of the reflection equation different from 
the factorizability requirement is related with the covariance properties of the 
X-matrix: if K satisfies the RE then the same is true for the transformed 
K' = TKT" 1 . It follows from (E3), (28) 



R X2 (u, v)YK 1 {u)TR 2l (v, -u)K 2 (v) = (33) 

= K 2 (v)R 12 {u, -v)VK x (u)TR 21 (-v, -u). 
Further generalizations of the YBE and the corresponding RE [37] are 



related to the colored algebras |59] and/or more complicated commutation 
relations among entries of the T-matrices and the i^-matrices. However the 
structure of these equations can be easily obtained following the standard ZF 
algebra pattern. In particular, if the multiplicative factor u(a, b) of the color 
ZF algebra is nondegenerate ||59||, where a, b G A and A is an abelian grading 



group, then the graded FRT-relation is 

R 12 (u, i;)r- 1 T 1 (u)TT 2 (v) = T 2 (v^T, (u)TR 12 (u, v) (34) 
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where (T) ab . cd = co(a, b)S ac 5bd and we identify for simplicity the matrix indices 
with the grading group elements (which is not always the the Z 2 - 

grading). Now T is not necessary unipotent. The corresponding RE follows 
from the covariance arguments 

R 12 (u, v)Y- 1 K x (u)TR 21 (v, -u)K 2 (v) = (35) 



K 2 (v)R 12 (u, -v)!- 1 ^ {u)TR 2l 



-u . 



To extract the commuting functionals of the ^4-graded monodromy matrix T 
or the corresponding Sklyanin's matrix T with the entries as the homogeneous 



elements of the grading the ^4-graded trace |59| has to be used 

t(u) = tr A T(u) = ^u(a,a)T aa (u) . 
This trace can be considered as a particular example of the quantum trace pl| 



with T as the i?-matrix. It is also possible to extend these constructions further 
by using instead of the multiplicative factor u(a,b) and/or T an appropriate 
i?-matrix, however the main problem for the moment is to find an interesting 
example the solution of which requires the mentioned above equations. 



5 Integrable models with finite degrees of free- 
dom. 

The intensive development of the YBE and the quantum group theory was 
strongly influenced by the conformal field theory. The latter one as most 
of the field theoretical integrable models has the two dimensional space-time. 
However, the rich structure of the YB-algebra and RE-algebra or the quadratic 
i?-matrix algebras with the spectral parameter dependence permits to include 
into this formalism variety of known integrable models with finite degrees of 
freedom (e.g. [38-40, 41, 46]) and to find new ones, which are physically 
interesting systems in the space of the three (and more) dimensions. 

Let us consider as an example the Kowalewski-Chaplygin-Goryachev top 
(KCG top). The Hamiltonian of this model is [46]. 



H = ^(J x 2 + J\ + 2J 3 2 ) + c 1 x 1 + c 2 x 2 + c 3 (x 2 - x 2 2 ) + c 4 xix 2 + c 5 /xl, (36) 

where arbitrary constants and Jj, Xi are the angular mo- 

menta and coordinates. This system is integrable provided the constraint: 

3 

l = J2xjJj = 0- (37) 
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One gets the famous Kowalewski's top for c 3 = C4 = C5 = 0. The corre- 
sponding auxiliary linear problem (the L-operator) is related to the simplest 
i?-matrix: the Yang solution (9) for sl(2) with c = k. The L-operator has 
quadratic dependence on the spectral parameter u 



where yf = y± ±y 5 and the same for yf; y , y 7 are eight dynamical variables 
of the model. 

The corresponding YB-algebra, generated by four entries a(u), d(u) of 
the L-operator, has its centre generated by the g-determinant of the L-operator 
©[5,3]: 



det q L(u) = a(u + in/2)d(u — in/2) — b(u + in/2)c{u — in/2) . . 

= Qiu 2 + Q 2 u — Q 3 + k 2 /AQi , {6J} 

where Q = y , Qi = y yz - y\ - yh Q2 = 2/22/3 - 2y 4 2/e - 2y 5 y 7 , Q 3 = 
V6 + V7- |{2/i> 2/2} + l^yoVs, {Vi, Vj} = ViVj + VjVi- 

The YB-algebra for the entries of T(u) = L(u) results in quadratic com- 
mutation relations for the dynamical variables y^. The nontrivial problem is 
to realize the latter ones in terms of physically significant variables. One of 
the realization is given [7] by the momenta Pi and coordinates q^, % — 1,2, of 
the two site Toda lattice, so that the L-operator (|3"8D is the product of two 
elementary ones (i = 1, 2) 

Z *M=(e-» ft o" 6 ") • (40) 

The realization we are looking for defined by the generators Jk, Xk, k = 
1, 2, 3 of the Lie algebra e(3), provided that I = J2k=i x kJk = [46] 



2/o = l y A = ibxi 

yi = -(Jf + J| + i + 2a/xf) y 5 = ibx 2 

2/2 = -2L 3 2/6 = -hib{x 3 , Li} 

2/3 = b 2 xl y 7 = -±ib{x 3 , J 2 } 



(41) 



where a and b are constants related to the YB-algebra central elements. The 
integrals of motion are generated by the trace of (|38|) 

t(u) = trL(u) = a(u) + d(u) = u 2 - 2J 3 u — 2H — ^ , (42) 

where H = |(L^ + J| — b 2 x 2 ) + ajx\ is the Hamiltonian of the Newmann's 
systems. The solution of this system by the separation of variables (from 
Jk, Xk to the new ones) in the framework of the QISM is achieved using the 
functional Bethe Ansatz [7, 41]: the introduction of the new variables as the 
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operator roots Uj, i — 1, 2, of the entry c(-u) (^) and the conjugated variables 
rhi as the values of a{u) and <i(u) at these roots. Then, the eigenvalue equation 
for the transfer matrix t(u) (two degrees of freedom for the Newmann's case) 
is reduced to the one-dimensional problems. To embed the KCG top into such 
approach one has to use the RE and the Sklyanin monodromy matrix. 

The c-number solutions to the RE with the s/(2) i?-matrix have the form 

K a (u) = a a a + u(P~a + + /?+cr_), (43) 

where a = ±, (/3l, fit, = (1, —fti, P2, — 1)- Using the L-operator (RSTj 



one can construct according to the general recipe [29] the monodromy matrix 

T(u) = L(u)K4u - tx/2)a 2 L T (-u)a 2 (44) 
and the correspondent generating function t(u) of the integrals of motion 



t(u) =trK+(u + ix/2)T(u) . (45) 

The latter one gives rise to the KCG top Hamiltonian fl5Ep. 

Another wide class of models embedded recently into this scheme [49] in- 
cludes the so called quasi- solvable models; some of them are physically relevant. 
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